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ABSTRACT 

The distribution of luminous radio galaxies in galaxy clusters has been observed to be 
concentrated in the inner region. We consider the role of dynamical friction of massive 
galaxies (M ^ 10^^'^ Mq), assuming them to be hosts of luminous radio galaxies, and 
show that beginning with a Navarro-Frenk- White density profile of a cluster of mass 
Mci ~ 10^^ M0 of concentration c ~ 5 and collapsing at z 1, the density profile of 
radio galaxies evolve to a profile of concentration c ~ 25, as observed, in a time scale 
of t - 3-5 Gyr. 
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1 INTRODUCTION 

Recent observations have shown that powerful radio sources 
in galaxy clusters are concentrated toward the cluster cen- 
ter. In a survey of 30 clusters of galaxies, Morrison & Owen 
(2003) found that the spatial distribution of high luminosity 
radio galaxies in rich clusters is described by a small core ra- 
dius. In particular, they found that the spatial distribution 
of galaxies depend on their radio power — high luminosity 
radio galaxies have a smaller core radius (~ 0.12±0.02 Mpc) 
than low luminosity ones (~ 0.4 ± 0.08 Mpc). Recently, Lin 
& Mohr (2007) have collected data for 188 rich clusters, and 
they have analyzed the spatial distribution of radio galaxies 
contained within them. They found that the spatial distri- 
bution of powerful radio sources could be fit by the 'univer- 
sal' Navarro-Frenk- White (NFW) profile (Navarro, Frenk, 
White 1997) with a larger concentration parameter than 
needed for fitting the total mass distribution of the cluster. 
And similar to the findings to Morrison & Owen (2003), they 
discovered that the concentration of galaxies in the spatial 
distribution increased with their radio power. 

Many explanations can be put forward for such a steep 
spatial distribution of radio galaxies in clusters. It is possible 
that the increased number density of galaxies in the inner 
region facilitates mergers and interactions that can trigger 
active galactic phenomena and increase the radio luminosity 
of galaxies. Also, radio galaxies are mostly associated with 
giant ellipticals which generally inhabit the central regions 
of rich clusters. Moreover, the increased pressure of the hot 
X-ray gas of rich clusters in the inner region may enhance 
the radio luminosity of FRII radio galaxies by confining the 
lobes with high pressure. 

There is also a dynamical aspect of the problem besides 
these factors. If the hosts of radio loud objects in general are 
massive galaxies, then dynamical friction will lead to mass 
segregation in galaxy clusters, and radio galaxies will settle 



toward the cluster center over many orbits. The situation is 
analogous to the phenomenon of large density of millisecond 
pulsars in the central regions of globular clusters. This con- 
centration is believed to be due to favorable conditions in 
the cluster core for accretion spin-up of neutron stars lead- 
ing to the phenomenon of millisecond pulsars. However, it 
is first the process of mass segregation in globular clusters 
owing to dynamical friction that makes massive stars settle 
toward the center, thereafter leading to an increased num- 
ber density of neutron stars in the inner regions (Meylan & 
Heggie 1997). 

Mass segregation by the process of dynamical friction 
is also expected to occur in galaxy clusters and it will lead 
to a preponderance of massive galaxies toward the cluster 
center. This process is likely to occur simultaneously with 
other processes of galactic evolution that may make galaxies 
in the inner region prone to becoming a radio loud source. 
To study the dynamical aspect of the evolution of spatial 
distribution of radio sources, one would need to isolate it 
from other evolutionary process and determine if it is a sig- 
nificant process by itself. We wish to address this issue in 
this paper. 

The main motivation to isolate the dynamical aspect 
of the problem comes from the fact that the hosts of ra- 
dio galaxies are often massive galaxies, irrespective of their 
membership in clusters or groups of galaxies. Recent sta- 
tistical studies of radio loud objects show that their hosts 
are predominantly massive objects (containing black holes of 
mass ~ 10*-^ Mq) (Best et al. 2005). It is therefore reason- 
able to expect that dynamical evolution of massive galaxies 
in clusters will have a significant eff'ect on the spatial distri- 
bution of radio loud sources in galaxy clusters, irrespective 
of other processes that contribute toward increasing the ra- 
dio output of a galaxy. 

Here we focus on the dynamical evolution of the distri- 
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bution of massive galaxies in clusters from dynamical fric- 
tion, and compare with recent observations. 



2 HOSTS OF RADIO GALAXIES 

Powerful radio sources (with Piaghz ^ 10^** W Hz~^) have 
been found to be associated with giant ellipticals, the oldest 
and the most massive galaxies in the universe. Best et al. 
(2005) have discussed the mass dependence of radio activ- 
ity with the help of a statistical analysis of the radio prop- 
erties of galaxies from the 2dF survey. They have shown 
shown that the probability of a galaxy containing a central 
black hole of a certain mass to have a given (or larger) ra- 
dio luminosity depends strongly on the black hole mass; the 
probability scales as oc Mg%. The (volume weighted) dis- 
tribution of the black hole mass in the sample showed that 
most radio loud objects had a black hole mass in the range 
Mbh ~ 10**-10^ Mq, peaking at Mbh ~ 10**'^ Mq. The 
total mass of the host galaxy of radio loud sources is related 
to the central black hole mass. Fine et al. (2006) found from 
the analysis of the 2dF QSO survey that for a QSO host 
containing a central black hole of mass Mbh ~ lO**-* Mq, 
the total dark matter mass of the host galaxy is ~ 10^^'^ 

Mq. 

Based on these considerations, we will consider galaxies 
of total (dark matter) mass of ~ 10^^'^ Mq, and containing 
a central black hole of mass Mhh ~ 10*'^ Mg, as examples of 
host galaxies of radio loud objects in our calculation below. 



3 DYNAMICAL FRICTION 

Consider a galaxy cluster of total mass Md and velocity 
dispersion a. We assume that the average mass distribution 
inside the cluster follows the Navarro-Frenk- White (NFW) 
profile (Navarro, Frenk, White 1997), p = ps/ {{r /rs)(l + 
(r/ra))), where = Vyir/c, c is the concentration pa- 
rameter, and Tvir is the virial radius. The virial radius is 
fixed by the overdensity estimated from spherical collapse 
model, which is approximately A{z = 0) ~ 100 at the 
present epoch, for the standard ACDM cosmological model 
(Om = 0.3, r^A = 1 — f^m) (see, e.g., Komatsu & Seljak ). 
The characteristic density ps is then given by. 



A(z) pc{z) 



3[ln(l + c) 



(1) 



where Pc{z) is the critical density of the universe at redshift 
z. For a large range of length scales (barring the central and 
outer regions of the cluster), the mass density p cx r~^. The 
mass distribution in these parts will be characterized by a 
near constant velocity dispersion a. 

We assume a concentration parameter c — 5 for a rich 
cluster as suggested by N-body simulations (Navarro et al. 
1996). We also assume a Maxwellian distribution of veloci- 
ties of objects in the clusters with dispersion a, and assume 
the average mass of galaxies to be m. 

Consider the motion of a massive galaxy of mass M and 
speed V in this cluster. Its motion in the background poten- 
tial of numerous small galaxies of mass m will be perturbed 



by random interactions with these galaxies. These pertur- 
bations will amount to a force of dynamical friction given 
by, 



drag 



erf(X)- 



2X 
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where n is the number density of galaxies of mass m, and 
X — V/^/2a. The dynamical friction will slowly change its 
orbital parameters and the massive galaxy will slowly settle 
toward the cluster center over many orbits. 

Consider a near circular orbit of the massive galaxy. 
Then, for a large range of length scales in the cluster, we 
have V = \/2cr, and X = 1. The coulomb logarithm is ap- 
proximately In A ~ In '''"gXf^ (Binney & Tremaine 1987, 
eqn 7-13b). For M ~ 10^^'^ Mq, and V ~ 1300 km/s (see 
below for the choice of this value), we have In A ~ In yf^^. 
If we take bmax as the cluster core radius which is ~ 100 kpc 
for rich clusters, then In A ~ 2.5. In this case, the product 



-X' 



In A ~ 1.1, and we assume it to be 



[erf(X)-2|exp(- 

unity for simplicity. Also, we will use the average density of 
matter p = nm in the above expression. 

If we consider an inspiraling orbit of a massive galaxy 
owing to the dynamical friction from its interactions with 
average mass galaxies, then one can estimate the time scale 
over which the distance of the massive galaxy from the 
cluster center would significantly decrease. Since the angu- 
lar momentum L ~ MVr changes according to dL/dt — 
Fdrag X r, wc havc 



(3) 
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since the circular velocity V changes little in the region of the 
cluster potential under consideration (where p oc and 
the circular velocity V = \/2o")- Ostriker & Turner (1979) 
also used this equation to study dynamical friction in galaxy 
clusters (see also Nusser & Sheth 1999). 

We can estimate the time scale of dynamical friction as. 



Mr. 



Uy„ = r/f ~ j^l^. For V ~ 10^ km/s, M ~ 10^ iviq, 
and using the matter density at the characteristic radius rs 
(using eqn[l}, ps ~ iq-^^-^ g cm~'^, for c = 5 and z — 0, we 
have. 
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V 1012 5 Mq^ V 10-26-6 g/cc^ 



(4) 



It is believed that rich clusters formed at 2 ~ 1, judging 
from the lack of evolution in X-ray luminosity function of 
clusters up to 2; ~ 1 (see, e.g., Rosati et al. 2002). The 
mass distribution would have settled into a NFW profile 
by that epoch. If such a cluster contained a few massive 
galaxies far from its center, then the above estimate shows 
that these massive galaxies would have spiralled toward the 
cluster center in a few Gyr time scale. 

The time evolution of the orbit of a massive galaxy em- 
bedded in a galaxy cluster can be obtained from numerically 
solving eqn|3] A few representative cases are shown in Figure 
1 in for M = lO^^ '^ Mq,V = 10^ km s"^ We compute the 
density distribution of the cluster halo for M^j = lO'^^ Mq 
collapsing at z = 1. We emphasize that the results shown 
in Figure 1 are approximate and are based on several as- 
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Figure 1. The evolution of orbit for a massive galaxy (M = 
3 X 10-"^^) embedded in a rich cluster (Md = 10^^ M0) is shown, 
for a few initial radii. 

sumptions. Firstly, we have assumed circular orbit. But the 
results shown in Figure 1 can be interpreted as the evolu- 
tion of the radial distance averaged over an orbital period. 
Also, we have neglected the effect of the massive galaxy be- 
ing an extended object. In reality, the galaxy will be tidally 
stripped, decreasing its mass, as it spirals inward (see Nusser 
and Sheth (1999) for a discussion on this effect). 

We can then ask how the distribution of such massive 
galaxies would evolve in a statistical sense, beginning with 
the NFW profile of a given concentration parameter. A given 
cluster will contain only a handful of massive galaxies, but 
we can consider a statistical ensemble of clusters with the 
same concentration parameter and total mass, which will 
contain several massive galaxies with initial orbital param- 
eters such that the total mass distribution of the ensemble 
cluster resembles a NFW profile of the assumed concentra- 
tion parameter (here c = 5). Slowly, over time and many 
orbits, the massive galaxies will settle toward the center, 
and their mass distribution will change. 

To estimate the gradual change in the mass distribu- 
tion of massive galaxies immersed in a gravitational po- 
tential determined by several average mass galaxies, we di- 
vide the cluster into several annuli, and compute the cluster 
halo density in the i-th annulus as that given by the NFW 
profile for the average radius in that annulus. Then we as- 
cribe an initial number of massive galaxies in each annulus, 
NM{i,t — 0), so that the density profile of these galaxies 
{NM(i,t — 0)/vol{i), where vol{i) is the volume of the i-th 
annulus) follows the NFW profile with same concentration 
(c = 5 here). We normalize the number of massive galaxies 
in different annuli by the density at the outermost annulus, 
which does not change owing to negligible dynamical friction 
at that radius. 

We first calculate the change in radius of these massive 
galaxies initially stationed at each annulus using eqn (|3]l . We 



compute their radii after an elapsed time interval, and then 
we re-calculate the number of massive galaxies in all annuli 
NM{i,t) at this epoch, again normalized by the number of 
massive galaxies in the outermost annulus. From this we 
calculate the number density of massive galaxies in each 
annulus as a function of time. The result will depend on 
the mass assumed for the massive galaxies and the time 
elapsed. Our previous estimates show that, we can expect a 
significant evolution in the distribution of massive galaxies 
over a time scale of a few Gyr for galaxies with M ~ 10^^'^ 

We consider a cluster of mass M^i = 10^^ M0 and con- 
centration parameter c = 5. The virial radius for a cluster 
collapsing at 2 ~ 1 is calculated to be r^ir ~ 2.6 Mpc. The 
appropriate velocity dispersion for such a cluster is taken 
from the scaling observed in simulations by Evrard et al. 
(2008) to be cr ~ 9.5 x 10* cm s"^ (for Hubble constant 
H = 70 km/s/Mpc); we assume it to be independent of lo- 
cation inside the cluster. We consider massive galaxies with 
M = lO^^ '^ Mq and V = \f2a ~ 1.3 x 10* cm s'^ 



4 RESULTS 

We show in Figure 2 the initial mass distribution by dotted 
line; it is a NFW mass profile with c = 5, for Md = 10^^ 
M0 collapsing at z = 1. We plot the radius as a fraction of 
r2oo- Then, the mass distribution of an ensemble of galaxies 
with Af = 3 X 10^^ Mq is shown after t = 1.5, 3 and 5 Gyr 
with solid lines. It is found that the distribution steepens 
with time. The distribution at outer radii does change much, 
but massive galaxies in the middle region settles toward the 
center, increasing the density of these galaxies in the bins of 
radii in these regions. The dashed line shows a NFW profile 
with c = 25. All profiles have been drawn keeping the mass 
density at the outermost bin fixed in time, for comparison. 

We therefore find that the distribution of massive galax- 
ies in a statistical ensemble of clusters steepens with time. 
The steepening depends on the assumed mass of the mas- 
sive galaxies and the time elapsed. In our fiducial case of 
M = 3 X 10^^ M© immersed in a cluster of total mass 
Mci -\r 10^^ Mq collapsing at z = 1, we find that a dis- 
tribution with initial c = 5 has changed to a distribution 
with c = 25 over a time scale t ~ 3-5 Gyr. 

It is interesting to compare this time scale with the 
cosmological lookback time. The look back time to redshift 
z ~ 0.5 is t ~ 4 Gyr. If rich clusters formed at these red- 
shifts, then their mass distribution would have settled into 
the NFW profile by that epoch. If such a cluster contained 
a few massive galaxies far from its center, then it is conceiv- 
able that these massive galaxies would have spiralled toward 
the cluster center by the present epoch. 



5 DISCUSSIONS 

Lin & Mohr (2007) have stacked the data from 188 clusters 
containing 16,646 radio loud objects (with radio luminos- 
ity P > 10^"^ W Hz~^) to produce the density profile of 
these objects. The number of radio loud objects within the 
radius r2oo in their study was 836. They found that the den- 
sity profile corresponds to a NFW profile with c ~ 25 ± 7, 
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Figure 2. The density profile of galaxies is plotted here. The 

dotted line shows the initial profile, and the solid lines correspond 
to the density profile of massive galaxies with M = 3 X lO-*^^ Mq 
after t = 1.5, 3 and 5 Gyr, the steepest ijrofilc being for the longest 
time interval. The dashed line corresponds to a NFW profile with 
c = 25 (normalized to have the same density as other profiles at 
the outermost bin). 

excluding the brightest cluster galaxy (BCG). With the in- 
clusion of BCGs which they defined as radio loud objects 
within 0.1r2oo, the profile steepens to c ~ 521^4- 

The solid curve in Figure 2 steepens further below a ra- 
dius 0.1r2()(), and excluding this part (equivalent to excluding 
the statistics of BCGs), we find that the density distribu- 
tion is comparable to a NFW profile with c ~ 25, close to 
the observed value. It then appears the distribution of radio 
loud objects in the ensemble cluster (produced by stacking 
the data of numerous clusters) can bo explained by assum- 
ing that the original mass distribution was a standard NFW 
profile with c ~ 5, and that massive galaxies slowly sink 
toward the cluster center owing to dynamical friction. The 
time taken to match the observed distribution ist ^ 3-5 Gyr 
if the massive galaxies have M ~ 10^^'® M©, and belonging 
to a cluster of Md = W^^ M©. 

It is interesting to note that Lin & Mohr (2007) have 
found that more powerful radio sources are more centrally 
concentrated than the weaker ones. For example, they found 
that radio sources with luminosity larger than 10'^^"'' W 
Hz~^ have a distribution comparable to a NFW profile with 
c ~ 59 ± 11. This trend is expected in the above scenario if 
the probability of a galaxy to have a given radio luminosity 
increases with its total mass, which indeed seems to be the 
case. Best et al. (2005) found from the 2dF survey data that 
the probability of a galaxy with a black hole mass Mbh to 
have radio power larger than a given value approximately 
scales as M^fj . If the central black hole mass scales linearly 
with the dynamical mass May„ of the host galaxy as (e.g., 
Hopkins et al. 2007), then one expects the hosts of more 
powerful radio sources to be more massive galaxies. 



We note here that dynamical friction cannot be the 
complete story behind the concentration of radio loud ob- 
jects in clusters. Masses of galaxies alone do not determine 
their radio properties; there are many factors that contribute 
to its radio luminosity, like mergers and interactions with 
other galaxies, and properties of the medium confining ra- 
dio lobes. Here we have isolated the dynamical aspect of the 
issue, and the results show that dynamical friction can also 
be as important as other factors in producing a steep density 
profile of radio sources in clusters. 



6 SUMMARY 

We study the dynamical friction of massive galaxies — which 
can often be hosts of radio loud objects — in galaxy clusters, 
and find the the time scale of mass segregation to be t ^ 3 5 
Gyr, for galaxies of mass ~ 10^^'^ M© within clusters of 
mass Mci ~ 10^^ M©. We show that this effect can help 
explain the observed distribution of radio loud objects in 
rich clusters of galaxies. 
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